Evolution of topological order in Xe films on a quasicrystal surface 
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We report results of the first computer simulation studies of a physically adsorbed gas on a quasicrys- 
talline surface, Xe on decagonal Al-Ni-Co. The grand canonical Monte Carlo method is employed, 
using a semi-empirical gas-surface interaction, based on conventional combining rules, and the usual 
Lennard- Jones Xe-Xe interaction. The resulting adsorption isotherms and calculated structures are 
consistent with the results of LEED experimental data. The evolution of the bulk film begins in 
the second layer, while the low coverage behavior is epitaxial. This transition from 5-fold to 6-fold 
ordering is temperature dependent, occurring earlier (at lower coverage) for the higher temperatures. 



The observed unusual electronic [1,2] and frictional 
[3-5] properties of quasicrystal surfaces stimulate inter- 
esting fundamental questions about how these and other 
physical properties are altered by quasiperiodicity. Re- 
cent progress in the characterization and preparation of 
quasicrystal surfaces raises new possibilities for their use 
as substrates in the growth of films having novel struc- 
tural, electronic, dynamic and mechanical properties [6], 
or as templates to facilitate the formation of well-ordered 
and/or size-selected nanoscale aperiodic arrays of clusters 
or molecules [7,8]. 

The physical behavior of systems involving competing 
interactions in adsorption is a subject of continuing in- 
terest [9] and is particularly relevant to the growth of 
thin films. Several different growth modes have been ob- 
served for the growth of metal films on quasicrystals. For 
instance, Cu on Al-Pd-Mn grows layer-by-layer to pro- 
duce a film having aperiodic order in one direction [10]. 
Al on Al-Cu-Fe forms small size-selected nanoclusters at 
very low coverages [11]. Ag on Al-Pd-Mn nucleates at 
strong-binding sites and ultimately grows as pyramidal 
hexagonal nanocrystals [12]. Sb and Bi on Al-Pd-Mn 
and Al-Ni-Co form quasiperiodic monolayers [13]. Some 
of these structures convert to different, often periodic, 
structures when annealed, and the issue of true equilib- 
rium in the film is sometimes complicated by a tendency 
for chemical intermixing with the substrate. The wide 
range of behavior observed so far indicates that, even in 
the absence of intermixing, film growth is strongly af- 
fected by chemical interactions between adsorbate and 
substrate. In order to separate these chemical effects 
from those specific to quasiperiodic order, we have stud- 
ied the adsorption of rare gases on a quasicrystal surface, 
where both the gas-gas and gas-surface interactions are 
believed to be simple, i.e., appreciable chemical interac- 
tions and adsorbate-induced surface reconstructions are 
absent. 

In the present work, we explore the implications of 
structural mismatch by evaluating the nature of Xe ad- 



sorption on a quasicrystal substrate, namely the 10-fold 
surface of decagonal Al-Ni-Co. This study of thermal 
and structural properties employs grand canonical Monte 
Carlo (GCMC) simulations, with which we have extensive 
experience [14-17]. The calculations employ the same 
potential function we used earlier to compute the low 
coverage adsorption with the virial expansion [18]. Us- 
ing the GCMC method, we compute the film properties 
for specified thermodynamic conditions. A hard wall at 
10 nm above the surface is used to confine the coexist- 
ing vapor phase. We take a square section of the surface 
A, of side 5.12 nm, to be the unit cell in the simulation, 
for which we assume periodic boundary conditions. This 
approach sacrifices accuracy of the long range QC struc- 
ture. However, such a simplification is numerically useful 
for these simulations. Since the cell is large relative to 
the Xe size it is accurately representative of order on 
short-to-moderate length scales. The simulation results, 
presented below, are remarkably consistent with both the 
results from the virial calculations [18] and with our ex- 
periment [19]: a monolayer film is found, the ordering 
of which reflects many aspects of the underlying qua- 
sicrystalline structure. During the onset of second layer 
adsorption, there is a transition into a 6-fold structure, 
which after several layers is consistent with a bulk Xe 
structure. 

The Xe-surface potential used here, shown in Figure 1, 
was calculated earlier by first developing an empirical ad- 
sorption potential, based on simple combining rules and 
known properties of bulk phases of solid Al, Ni and Co 
crystals [18]. The potential is based on a summation of 
two-body interactions between the Xe and the individ- 
ual constituent atoms of the substrate: Al, Ni and Co. 
The gas-gas potential is taken simply to be a Lennard- 
Jones (LJ) interaction, with parameter values for Xe: 
(7 = 221 K and e = 0.41 nm. The Xe-substrate pair 
interactions are also assumed to have LJ form, with pa- 
rameter values taken from traditional combining rules, 
using atomic sizes derived from bulk crystalline lattice 
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constants [18-20]. In the calculation of the adsorption 
potential, we assumed a structure of the virgin surface 
taken from an empirical fit to LEED data [211. 
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FIG. 1. Computed potential energy for Xe on AINiCo 
and potential energy scale (right), obtained by minimizing 
V(z,y,z) with respect to z variation. 

The Xe adsorption potential derived with this proce- 
dure is both deep and highly corrugated. Depending on 
the lateral position (x, y) across the surface, the max- 
imum depth as a function of normal coordinate (z) is 
typically in the range D(x,y) ~ 150 to 250 meV. This 
potential is called "deep" because the record maximum 
well-depth for Xe on a periodic surface is about 160 meV, 
viz. on graphite [22]; the record minimum well-depth is 
about 28 meV, on Cs [23]. While the term "corruga- 
tion" is not well-defined for a non-periodic surface, one 
estimate of its magnitude comes from the standard devi- 
ation of the laterally varying well-depth, which is about 
50 meV on this surface. Thus the well-depth's fluctu- 
ation is about 25% of its average magnitude, which is 
sufficiently large as to warrant the description "highly 
corrugated", in our opinion. Using this potential, we 
computed the adsorption properties at low vapor pres- 
sure, P, using the virial expansion (including the first 
three terms). The results of that analysis were found to 
be semi-quantitatively consistent with our experimental 
data in the low coverage regime [18]. 

Figure 2 shows adsorption isotherms computed with 
the GCMC simulations. The temperature range explored 
here extends from 70 K to 286 K (the triple temperature 
of Xe is 161.4 K). The plotted quantity is the thermody- 
namic excess coverage, N x /A, defined as the difference 
between the total number of atoms in the simulation 
cell, and the number that would be present if the cell 
were filled with uniform vapor at the specified values of 
P and T. Detailed inspection of the isotherms reveals 
that there is continuous film growth (i.e. complete wet- 
ting) for all temperatures above the triple point. This 
behavior persists to some temperature below the triple 
temperature but we have not clearly established where 



the wetting transition occurs. At 77 K, at least 5 layers 
form before the onset of bulk condensation. The inset in 
Figure 2 shows the density profile, p(z), in the direction 
perpendicular to the surface, for a total surface coverage 
of 26.28 atoms/nm 2 , corresponding to point "e" on the 
77 K isotherm. One observes that the (perpendicular) 
layer structure is well defined, in spite of the large lateral 
variation of the potential. 
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FIG. 2. Computed isotherms from 70 K to 280 K in steps of 
10 K. Two additional isobars at 77 K and 286 K are shown. 
The isotherm with solid circles is for 77 K, while the high- 
est-temperature isotherm is for 286 K. The inset shows the 
density profile p(z) for P = 0.259 x 10" 4 atm and 77 K, indi- 
cated by point (e). 

Figure 3 shows the density variations within the top 
layer for several points on the 77 K isotherm as speci- 
fied in Figure 2. At the lowest pressure (point a), one 
observes a film that has atoms localized in the deepest 
parts of the potential, and the Fourier transform (FT) 
of that density function is 10-fold symmetric, reflecting 
the substrate symmetry. At the top of the first-layer step 
of the isotherm (point b), the density variation is more 
uniform, now having points of localized 5-fold and 6-fold 
symmetry. The FT of this structure is still 10-fold, i.e. 
still reflecting the substrate symmetry. By the time the 
adsorption proceeds to point c, however, the layer has 
more points of local 6-fold symmetry, and the FT dis- 
plays 6-fold symmetry. When the full bilayer has formed 
(point d) the in-plane symmetry is clearly 6-fold locally, 
with some dislocation defects, and the FT indicates an 
overall 6-fold symmetry of the structure. When the 6-fold 
structure forms, it is aligned so that the close-packed di- 
rection of the Xe is parallel to a principal 5-fold direction 
of the substrate surface. For any given simulation run, 
the selection of 5-fold/6-fold alignment is arbitrary, and 
in Figure 3d, it can be seen that there are actually two 
alignments present. 

Xe adsorption on this surface was studied earlier using 
low-energy electron diffraction, and isobar measurements 
indicate that the Xe film grows layer- by- layer in the tem- 
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perature range 65K to 80K [19], consistent with the sim- 
ulations described above. Figure 3 shows the LEED pat- 
terns obtained from adsorbed Xe on decagonal Al-Ni-Co 
under similar conditions to the simulations. At the low- 
est coverage, the only discernible change in the LEED 
pattern from that of the clean surface is an attenuation 
of the substrate beams. After the adsorption of one layer 
(coverage determined by the isobar measurements [19]) 
there are still no resolvable features that would indicate 
an over layer having order different from the substrate. 
At the onset of the adsorption of the second layer, how- 
ever, the LEED pattern shows new diffraction spots that 
correspond to 5 rotational domains of a hexagonal struc- 
ture. Within each of these domains, the close-packed 
direction of the Xe is aligned with the 5-fold directions 
of the substrate [19], as observed in the simulation. In 
the experiments, all possible alignments are observed ow- 
ing to the presence of all possible rotational alignments 
present within the width of the electron beam (0.25 mm), 
and there is no evidence of step-pinning of the over layer. 
When the second layer is complete, these spots are well- 
defined and their widths are the same as the substrate 
spots, indicating a coherence length of at least 15 nm. A 
dynamical LEED analysis of the intensities indicates that 
the structure of the multilayer film is consistent with fee 
Xe(lll). 
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FIG. 3. Intralayer density plots for adsorption at 77 K, 
Fourier transforms of these plots and LEED patterns corre- 
sponding to similar conditions. The plots show the density in 
the top layer, corresponding to the points noted in Figure 2. 



While the experiments are restricted to a compara- 
tively narrow range in T and P (the LEED experiments 
require P < 10 _ 7 bar and feasible equilibration times 
require P > 10 -13 bar) the simulations are not so re- 
stricted. Therefore, we have also investigated the ad- 
sorption of Xe at a very low temperature (20 K) and a 
high temperature (160 K). The density plots and FT's 
at the coverage corresponding to point c in Figure 3 is 
shown on Figure 4 for the three temperatures, 20 K, 77 
K and 160 K. The trend observed in the ordering is that 
the onset of 6-fold ordering occurs earlier at the higher 
temperatures. Six-fold ordering is already present at this 
coverage for 160 K, but is barely present at 77 K and is 
clearly not present at 20 K. This trend is consistent with 
the substrate potential being a bigger influence on the 
film's structure at lower T, implying that the develop- 
ment of the bulk Xe structure is entropy- driven; this is 
a kind of wetting transition of the bulk Xe phase, at- 
tributable to a diminishing interfacial free energy cost 
with increasing T. More thermal disorder is also evident 
in the film structures at higher T. 
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FIG. 4. Density plots and Fourier transforms for monolayer 
films (corresponding to point c in Figure 2) at 20 K, 77 K, 
and 160 K. 

Interestingly, at all T studied, stacking faults are evi- 
dent in the multilayer films. Their origin appears to be 
dislocations in the layers, which are most prevalent at the 
highest temperatures studied, as expected for entropic 
reasons. This is consistent with x-ray diffraction studies 
of the growth of Xe on Ag(lll), where stacking faults 
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were observed for Xe growth under various growth con- 
ditions [24,25], although the overall structure observed 
was fcc(lll). Such a stacking fault is evident in Fig- 
ure 5, which shows a superposition of Xe layers 2 and 
4 at point e in Figure 2. The coincidence of the atom 
locations in the top left part of this figure is consistent 
with an hep structure (ABAB stacking) whereas the off- 
set observed in the lower left is consistent with an fee 
structure (ABC stacking). We note that while bulk Xe 
has an fee structure, and indeed an fee structure was 
found for the multilayer film in the LEED study, calcula- 
tions of the bulk structure using LJ pair potentials such 
as those employed here result in a more stable hep struc- 
ture [26] . This difference apparently arises from a neglect 
(using the pair potentials) of d-orbital overlap interac- 
tions, which are more effective in fee than in hep struc- 
tures [26]. This model agrees quantitatively with mea- 
surements of energy differences in fee and hep Ar near its 
triple point [27]. The presence of dislocations and discli- 
nations in the monolayer films and their variation with 
temperature raises interesting questions concerning the 
possibility of dislocation- mediated melting [28-30], but 
the limited simulation size precludes any such analysis. 
Larger simulations with Xe and different noble gases are 
intended for the future. 
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FIG. 5. Density plot corresponding to point c on Figure 
2, showing a superposition of the density slices for the 2 nd 
and 4 th layers. In the top right, 4^-layer atoms are located 
directly above the 2 nd layer atoms, whereas in other regions, 
such as lower left, the two layers are offset. 

In summary, our simulation results are consistent with 
the experimental data over the somewhat limited range 
explored thus far. A consistent pattern has emerged; the 
quasicrystalline order is transmitted to films of thickness 
less than about two layers, disappearing upon further 
growth of the film. The results of this study encourage 
confidence in the use of physical adsorption to probe film 
growth. The extension of these calculations to the study 
of other properties such as lattice dynamics and friction, 
as well as other adsorbates, would seem well justified. 
Such studies will yield a comprehensive understanding of 
competing interactions in physisorbed layers on aperiodic 
substrates, as found previously on periodic surfaces. 

This research was supported by NSF grant DMR- 
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